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ABSTRACT 

We study the properties of gas in and around 10^^ Mq haloes at ^ = 2 using a suite 
of high-resolution cosmological hydrodynamic ‘zoom’ simulations. We quantify the 
thermal and dynamical structure of these gaseous reservoirs in terms of their mean 
radial distributions and angular variability along different sight lines. With each halo 
simulated at three levels of increasing resolution, the highest reaching a baryon mass 
resolution of ^10,000 solar masses, we study the interaction of hlamentary inflow and 
the quasi-static hot halo atmosphere. We highlight the discrepancy between the spatial 
resolution available in the halo gas as opposed to within the galaxy itself, and hnd 
that stream morphologies become increasingly complex at higher resolution, with large 
coherent flows revealing density and temperature structure at progressively smaller 
scales. Moreover, multiple gas components co-exist at the same radius within the 
halo, making radially averaged analyses misleading. This is particularly true where the 
hot, quasi-static, high entropy halo atmosphere interacts with cold, rapidly inflowing, 
low entropy accretion. We investigate the process of gas virialization and identify 
different regimes for the heating of gas as it accretes from the intergalactic medium. 
Haloes at this mass have a weh-dehned virial shock, associated with a sharp jump 
in temperature and entropy at >1.25rvir- The presence, radius, and radial width 
of this boundary feature, however, vary not only from halo to halo, but also as a 
function of angular direction, covering roughly ^ 85% of the 47r sphere. Our hndings are 
relevant for the proper interpretation of observations pertaining to the circumgalactic 
medium, including evidence for large amounts of cold gas surrounding massive haloes 
at intermediate redshifts. 

Key words: methods: numerical - galaxies: evolution - galaxies: formation - galaxies: 
haloes - cosmology: theory 


1 INTRODUCTION 

Initially following the gravitational collapse of a dark matter 
overdensity, gaseous haloes subsequently grow through the 
accretion of baryons from the intergalactic medium (IGM). 
Their evolving structure across cosmic time has been the 
subject of theoretical as well as observational interest for 
several decades. As the transitional state between the diffuse 
IGM and the star-forming interstellar medium of galaxies, 
these gas reservoirs regulate the stellar growth of forming 


* E-mail: dnelson@cfa.harvard.edu 
t Hubble Fellow. 


galaxies. Understanding not only the structure of halo gas, 
but also its origin and subsequent evolution, is therefore 
essential for any comprehensive theory of galaxy formation. 

Although the accretion of gas will follow that of dark 
matter, the presence of additional physical processes includ¬ 
ing hydrodynamical forces and radiative cooling imply ad¬ 
ditional complexity for the acquisition of baryons. In the 
classic picture, gas accreting from the intergalactic medium 
will shock heat to the virial temperature of the halo. If 
the radiative cooling timescale is sufRciently long, it will 
then form a hot, pressure supported atmosphere in approx¬ 
imate equilibrium ( Rees Ostriker]|1977| Silk||1977| White 
Sz ReespQTS). Gooling can proceed, delivering gas into the 
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halo centre ( White fc Frenk|1991 1. For sufficiently low mass 
haloes, this timescale will be short enough that a ‘stable 
virial shock’ cannot develop, and gas accretion from the 
IGM will proceed as rapidly as dynamically allowed ( Birn- 
boim fc Dekel||2003a ). However, this theoretical foundation 
is largely based on a one-dimensional picture, while dark 
matter haloes and their gaseous counterparts are decidedly 
not spherically symmetric. Numerical hydrodynamical sim¬ 
ulations have indicated that galaxies can acquire their gas in 
a fundamentally different manner ( Katz et al.||2003 Abadi 


et al.|2003 Keres et aL|2005 Ocvirk et al.|20Q8 ). In partic¬ 

ular, finding that coherent, filamentary inflows can fuel star 
formation in a central galaxy while avoiding shock heating 
to the virial temperature. Such streams arise naturally from 
the topology of large-scale structure, particularly at high 
redshifts of z > 2 (Dekel et al. 2009 Agertz et al. 2009 


Keres et al.||2009 Danovich et al. 2012). 

Gas inflow has been studied in the context of many key 
questions in galaxy formation, most fundamental of which is 
perhaps its link to star formation ( Oppenheimer et al.|201(3 


Gabor fc Bournaud|2014 Sanchez Almeida et al.|20T4 ) and 

the process of quenching (Birnboim & Dekel 2003b| Birn- 
boim et ar]|2007| [Gabor fc Dave||2012| [Feldmann fc Mayer 


2015| Aragon-Galvo et al. 2014). The accretion of cosmo¬ 
logical gas will leave a fundamental imprint on the thermal 
and dynamical properties of the quasi-static halo gas, as will 


outflows from energetic feedback processes in galaxies (Put¬ 


man et al. 2012). By studying the properties of this halo 


gas we can investigate the interplay of inflows and outflows 
with metals (e.g. Shen et al.|2013 Hummels et al.|2013 Ford 
et al.||2014 ) as well as with neutral hydrogen (e.g. Faucher- 


Gigu^re et al.||2011| |Faucher-Giguere et al.||2014| IFumagalh 

et al.||2014|^ 

The state of gas in and around galaxy haloes has also re¬ 
ceived significant observational scrutiny in the local universe 
as well as at z~2, near the peak of the cosmic star forma¬ 
tion rate. At these high redshifts both hydrogen and metals 
are accessible as absorption signatures in sightlines towards 
background objects, enabling a probe of the interaction be¬ 
tween galaxies and the IGM in the vicinity at this period 
of rapid stellar growth. Efforts include the ‘quasars prob- 


ing quasars’ series (Hennawi et al. 2006 Prochaska et al. 

2013 

2014 

, the Keck Baryonic Structure Survey (Steidel 

et al. 

2010 

Rudie et al. 2012 2013 

Turner et al. 

2014b), 

and other 2 < 2 ; <3 studies (Simcoe et al.||2006 Fieri et al. 

2014 

Rubin et al.||2014 Grighton et al.||2015). Gollectively 


they probe the covering fractions, radial profiles, and kine¬ 
matics of HI and many metal ions including those of oxygen, 
carbon, neon, silicon, and magnesium. 

Perhaps the most puzzling discovery which remains at 
the present unreconciled with theory is the presence of a 
large amount of cold (~10^K), metal-enriched gas widely 
distributed at least out to the virial r adius of massive 
(10^^ ~ 10^^'^ M© ) haloes. For instance, Prochaska et al. 


(2014) find a covering fraction approaching unity for strong 


C II absorption at a hundred physical kiloparsecs or more 
from the halo centre. The origin of this gas and the process 
by which it is either maintained or replenished within the hot 
halo remains uncertain. Gonsequently, this regime provides 
a powerful test-bed for current hydrodynamic simulations 
of galaxy formation. Not only in terms of the physical mod¬ 
elling of feedback and the resultant galactic-scale outflows 


( Muratov et al.||2015 ), but more fundamentally in terms of 
our ability to resolve the gas-dynamical processes and the 
spatial scales relevant for the physics of cosmological gas 
accretion. 

This paper investigates the thermal and dynamical 
structure of halo gas in eight simulated ~ 10^^ M© haloes 
at z — 2. In Section we describe the simulation technique 
and analysis methodology. Section addresses the issue of 
resolution for halo gas and presents a visual overview of the 
systems. Section [^considers the radially averaged gas prop¬ 
erties, while Section expands this analysis to explore the 
angular variability of halo structure. We discuss our results 
in the context of observations of the gas content of haloes 
in Section [G] and summarize our conclusions in Section 0 


2 METHODS 


2.1 Initial Conditions 


All simulations evolve initial conditions which are a ran¬ 
dom realization of a WMAP-9 consistent cosmology (Qa,o = 
0.736, Qm,o = 0.264, 0^,0 = 0.0441, erg = 0.8 05, n, = o!967 
and h = 0.712). We use the MUSIC code (Hahn & Abel 
2011 vl.5, r375) to generate multi-mass ‘zoom’ IGs, under 
the 2LPT approximation and with a tabulated transfer func¬ 
tion from CAMB ( Lewis et al.||2000 ). First, we evolve a low 
resolution, dark matter only uniform periodic box of side- 
length 20/i“^ Mpe ~ 28.6 Mpe with 128^ particles (‘L7’, 
where LA'=2^), from a starting redshift of 2 : = 99 down to 
2 : = 2. At this redshift, there are 20 haloes with total mass be¬ 
tween lO^^'^M© and lO^^'^M© from which we choose eight 
at random to re-simulate at higher resolution. We do not 
select for any additional properties - e.g., merger history or 
environment. All particles within some factor of the virial 
radius of each selected halo (ranging from 3.6rvir to 7.0rvir 
in all cases, see Ohorbe et al. I poTIl i for relevant consider¬ 
ations) are identified at 2 ; = 2. We take rvir ='r 2 oo,crit the 
radius enclosing a mean overdensity 200 times the critical 
density. This factor was chosen by trial and error with eval¬ 
uation of contamination levels in low resolution test runs. 
The convex hull of the 2 ; = 99 positions of all selected par¬ 
ticles is then used to define the high resolution refinement 
region. 

For each halo, new initial conditions are generated for 
each of L9, LIO, and Lll, corresponding to 512^, 1024^, and 
2048^ total particles if the parent box were to be simulated 
at a uniform resolution. We note that the mass resolution of 
Lll is between Aquarius levels ‘3’ and ‘4’ (of e.g. Marinacci 
et al.|2014a Scannapieco et al.|2012 ), while the resolution of 


L9 is approximately equal to the resolution in modern, large 
volume cosmological simulations. Baryons are included by 
splitting each dark matter particle according to the cosmo¬ 
logical baryon fraction, into one DM particle and one gas 
cell, such that the centre of mass position and velocity are 
preserved. Therefore we do not consider a separate transfer 
function for the baryonic component. The files required to 
generate our initial conditions, including the CAMB transfer 
function, noise seeds, and convex hull point sets, are made 
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Table 1. General characteristics of our three resolution levels, L9, LIO, and Lll. First, the effective resolution of an equivalent uniform 
box. Next, the mean number of high resolution gas elements, number of timesteps, baryonic mass resolution, and dark matter mass 
resolution. The Plummer equivalent comoving gravitational softening lengths, and their physical values at z = 2. The minimum gas cell 
spatial size in physical parsecs at z = 2, and the mean gas cell spatial size in the halo, between O.Srvir and l.Orvir, in physical kiloparsecs 
at z = 2. 


Res 

Neff 

part 

NHR 

part 

At [#] 

^baryon [Mq ] 

iiidm [M© ] 

[pc] 

[pc] 

[pc] 

’■cen° [kpc] 

L9 

5123 

800,000 

80,000 

1.0 X 10 ® 

5.1 X 10® 

1430 

480 

31 

2.7 

LIO 

10243 

7,000,000 

260,000 

1.3 X 10 ® 

6.4 X 10® 

715 

240 

11 

1.6 

Lll 

20483 

64,000,000 

870,000 

1.6 X 10 ^ 

8.0 X 10"^ 

357 

120 

3.3 

0.8 


Table 2. Details on the eight simulated haloes: the total halo 
mass and (physical) virial radius at 2 ; = 2 from the parent box. 
The radius of the enclosing sphere at 2 ; = 2 used to define the 
Lagrangian region, in terms of Tvir, and the number of high reso¬ 
lution elements, for each of dark matter and gas. Both are listed 
for the Lll level only. The minimum radius from the halo centre 
reached by (i) contaminating low resolution dark matter parti¬ 
cles and (ii) Monte Carlo tracers originating in low resolution gas 
cells, in units of rvir- The total number of timesteps to reach 2 ; = 2. 


Halo 

halo 

^Par 
^ vir 

^Lll 

'hr 


_min 

'lr 

^min 

'lr 

At 

# 

[log M© ] 

[kpc] 

[rvir] 

[10®] 

[rvir] dm 

[rvir]tr 

# 

hO 

12.1 

114 

3.6 

70.0 

1.77 

2.16 

829714 

hi 

12.1 

104 

4.8 

66.7 

2.12 

2.75 

701681 

h2 

11.9 

92 

6.0 

24.2 

2.83 

3.02 

955189 

h3 

11.9 

96 

7.0 

33.9 

2.74 

3.23 

812983 

h4 

12.0 

103 

6.0 

68.4 

2.13 

2.89 

861224 

h5 

12.0 

103 

4.2 

59.9 

1.04 

1.04 

931088 

h6 

12.1 

97 

4.8 

74.4 

1.32 

1.59 

980918 

h7 

11.9 

94 

4.4 

52.8 

0.94 

1.93 

866242 


publicly available onlin^ The fundamental characteristics 
of each resolution level are given in Table while the phys¬ 
ical properties and numerical details for each of the eight 
haloes are detailed in Tabled 


2.2 Simulation Code and Physics 


We employ the AREPO code ( Springel|2010[ r25505) to solve 
the coupled equations of ideal continuum hydrodynamics 
and self-gravity. An unstructured, moving, Voronoi tessel¬ 
lation of the domain provides the spatial discretization for 
Godunov’s method with a directionally un-split MUSCL- 


Hancock scheme (van Leer 1977) and an exact Riemann 


solver, obtaining second order accuracy in space. Since we 
allow the mesh generating sites to move, with a velocity 
equal to the local fluid velocity field modulated by cor¬ 
rections required to maintain the regularity of the mesh, 
this numerical approach would be classified as an Arbitrary 
Lagrangian-Eulerian (ALE) scheme. Gravitational acceler¬ 
ations are computed using the Tree-PM approach, where 
long-range forces are calculated with a Eourier particle-mesh 


^ http://www.illustris-project.org/files/Nelsonl5b_ICs. 
zip 


method, medium-range forces with a hierarchical tree algo¬ 
rithm ( Barnes Hut|[l986 ), and short-range forces with 
direct summation (as in Springel||20Q5 ). A local, predictor- 
corrector type, hierarchical time stepping method obtains 
second order accuracy in time. Numerical parameters tan¬ 
gential to our current investigation - for example, related 
to mesh regularization or gravitational force accuracy - are 
detailed in Springel (2010) and Vogelsberger et al. (2012). 

We include a redshift-dependent, spatially uniform, ion¬ 
izing UV background field ( [Eaucher-Giguere et al. 2009). 
Gas loses internal energy from optically thin radiative cool¬ 
ing assuming a primordial H/He ratio (Katz et al. 1996). 
The production of metals and metal line cooling contribu¬ 
tions are not included.Star formation and the associated 
ISM pressurisation from unresolved supernovae are included 
with an effective equation of state modelling the ISM as a 
two-phase medium, following Springel & Hernquist (2003). 
Gas elements are stochastically converted into star parti¬ 
cles when the local gas density exceeds a threshold value 
of uh = 0.13 cm“^. Eurthermore, there is no resolved stel¬ 
lar feedback that would drive galactic-scale winds, nor any 
treatment of black holes or their associated feedback. All of 
the simulations considered in this work disregard the possi¬ 
ble effects of radiative transfer, magnetic fields, and cosmic 
rays. The set of implemented physics is essentially identical 
to the ‘moving mesh cosmology’ paper series as described in 


Nelson et al. (2013). 


We identify dark matter haloes and their gravitation¬ 
ally bound substructures using the SUBFIND algorithm 
(Springel et al. 2001 Dolag et al. 2009) which is applied 
on top of a friends-of-friends cluster identification. We refer 
to the most massive substructure in each FoF group as the 
halo itself. Merger trees are constructed using the SUBLINK 
code ( Rodriguez-Gomez et al.||2015 ) to link haloes to their 
progenitors and descendants at different points in time. 

All runs include the Monte Garlo tracer particle tech¬ 
nique (Genel et al. 2013) in order to follow the evolving 
properties of gas elements over time, with five tracers per 
initial gas cell. This is a probabilistic method, where tracers 
are exchanged between parents based explicitly on the cor¬ 
responding mass fluxes. By locating a subset of their unique 
IDs at each snapshot we can, by reference to their parents at 
that snapshot, reconstruct their spatial trajectory or ther¬ 
modynamic history. 


In the current work we make limited use of the assem¬ 
bly and accretion histories of the haloes, as determined with 
the tracer particles and merger trees, respectively. Specifi¬ 
cally, we verify the correspondence between haloes at differ- 
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Figure 1. The mass-weighted distribution of gas spatial resolu¬ 
tion as a function of radius, for all eight halos stacked together at 
Lll resolution. As a proxy for spatial resolution, we show the gas 
cell sphere-equivalent radii rceii- Each radial bin is normalized in¬ 
dependently, such that the colour mapping reaches its maximum 
intensity at each radius, independent of the radial distribution of 
gas mass. The dark vertical lines indicate l.Srvir, O.lSrvir and 
O.OTrvir, the last two bounding the radial range where a density 
bimodality arises due to the central galaxy. Substructures are ex¬ 
cised. 


ent resolution levels in order to comment on the convergence 
properties of our simulations. We reserve as future work - 
the second paper in this series - a quantitative analysis of 
the rates and modes of accretion and the impact of merging 
substructures. 


RESOLUTION CONSIDERATIONS AND 
VISUAL INSPECTION 


In the vast majority of galaxy formation simulations, spa¬ 
tial resolution is naturally adaptive and follows the hierar¬ 
chical clustering of structure formation in ACDM. This is 
true for the dark matter, where the Vlassov-Poisson equa¬ 
tions are solved with a Monte-Carlo approach. The N- 
body method is typically used (for notable exceptions see 
Yoshikawa et al.||2013| Hahn Angulo||2015 ). This natural 
adaptivity also holds for the gaseous component in particle 
methods like SPH, adaptive grid methods with the typi¬ 
cal density-based refinement criteria (for notable exceptions 
lapichino et al.||2QQ8 Rosdahl Blaiz^|2012 ), and the 


moving-mesh code used in this work. The resolution is there¬ 
fore better in collapsed structures than it is in under-dense 
regions, and maximal in galaxies themselves. The gaseous 
haloes surrounding galaxies are poorly resolved in compari¬ 
son, to the extent that modern cosmological simulations (e.g. 


Khandai et al.||2014| [Dubois et al.||2014| [Vogelsberger et ah] 


2014[ [Schaye et al.||2015 ) which just barely resolve the in¬ 
ternal structure of individual galaxies may be insufficiently 
capturing gas-dynamical processes in the halo. 

We begin by quantifying the spatial resolution of gas 
in the halo. Figure^ shows the mass-weighted two dimen¬ 
sional distribution of cell size rceii as a function of radius. 
As a proxy for the irregular shape of Voronoi polyhedra, 
we take the sphere-equivalent radius rceii = (3yceii/47r)^/^ 
where Well is the exact Voronoi cell volume. The mesh reg¬ 


ularization scheme ensures that this value is approximately 
equal to the actual geometrical distance between the gen¬ 
erating point and each face centre. All haloes are stacked 
together at the highest (Lll) resolution level, while black 
lines show the median relation for each halo separately. 

Beyond r/rvir > 1.25 the hydrodynamic resolution is 
nearly constant, but with a broad distribution spanning from 
~ 500 pc to ~3kpc (physical), depending on what level of 
overdensity the gas cell resides in. In the halo region, 0.25 < 
T’/rvir <1.0 the mean resolution scales from ~ 400 pc up to 
~ 1 kpc. The inner halo, typically bounded by r < 0.15rvir 
sees the mean cell size decrease rapidly under ~ 100 pc as 
dense gas features associated with the central galaxy begin 
to dominate. We return to the implications of the relatively 
coarse halo resolution in the discussion. 

In Figure we show a visual example of the resolution 
issue and a first look at the impact of higher hydrodynamical 
resolution in the halo regime. One halo (hO) is included at 
the three increasing zoom levels considered in this work, L9, 
LIO, and Lll, increasing from left to right. On the top row 
we plot a slice of the Voronoi mesh used to evolve the gas 
component 1^ The slice is centred on the central galaxy, which 
lies at the origin, and to focus on details we include only the 
upper-right quadrant. The outer red circle marks the virial 
radius of 114 physical kpc at z — 2. At all three resolution 
levels we see that the high-density gas near the origin is re¬ 
solved at the sub-kpc level and impossible to separate at this 
scale, while at the virial radius individual gas elements are 
substantially larger. At L9, characteristic of modern cosmo¬ 
logical simulations, the virial arc across a quadrant is only 
resolved by ~ 20 cells, while the radial structure between 
0.5rvir and l.Orvir is only resolved by ^10 cells. This is 
roughly consistent with a zeroth order estimate, assuming 
uniform density, for a number of cells per dimension of 


Acell,din 


m ) X Mhalo \ 

^cell J 


1/3 


50 


( 1 ) 


which are available to resolve the halo gas structure of a 
10^^ Mq halo. At Lll this increases to Aceii,dim — 250, but 
the scaling is obviously slow in three-dimensions. In the dis¬ 
cussion we consider a possible method for improving upon 
this resolution issue in the halo with future work. 

Across the bottom row we show a projection of mass- 
weighted gas temperature, where the colormap extends from 
cold, ~ 10^ K gas as blue to hot, ~ 10® K gas as red. Some 
interesting differences emerge as the gas resolution in the 
halo increases. The large cold inflow which begins to expe¬ 
rience substantial heating at (0.4 — 0.5) x rvir shows some 
temperature inhomogeneities at L9, but is essentially a sin¬ 
gle coherent structure. At LIO, this inflow is resolved into a 
number of smaller gas filaments by the time it crosses half 
the virial radius. At Lll the temperature structure becomes 
even more complex, throughout the halo, with small scale 
features emerging below the cell-size of the lower resolution 


^ Note that the intersection of a plane with the three- 
dimensional Voronoi tessellation does not in general produce a 
two-dimensional Voronoi graph. Therefore the shapes of individ¬ 
ual cells, as shown, do not reflect the actual regularity of the 
computational mesh. 
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Figure 2. Slice through the Voronoi mesh showing the underlying spatial discretization (top) and a projection through the linearly re¬ 
constructed, mass-weighted gas temperature field (bottom) with ~ 10^ K as blue and ~ 10® K as red. Halo hO is shown at 2 ; = 2. We include 
only the upper right quadrant, after centring the halo at the origin, in orthographic projection. The circles denote {0.15,0.5,1.0} rvir, the 
axes units are in physical kpc, and the projection depth is 238 kpc. The average gas cell size near the virial radius drops with increasing 
spatial resolution (from left to right), from ~2.7 kpc to ~0.8 kpc physical. Note that the exact state of gas in the halo, and the location 
of satellites, differs somewhat due to timing offsets between the three runs. 


runs. Clearly the morphology of accreting gas and its inter¬ 
action with the quasi-static hot halo material depends on 
numerical resolution to some degree. In the following sec¬ 
tions we further investigate this dependence by analysing 
the structure of halo gas at our highest available resolution. 


3.1 Visual Inspection 

We present a visual overview of the eight simulated haloes 
at z = 2, focusing on the virial scale of haloes 0-3 in Fig¬ 
ure and of haloes 4-7 in Figure one halo per row. In 
each case, we show mass-weighted projections of gas den¬ 
sity, temperature, entropy, and radial velocity in each col¬ 
umn (left to right). All gas cells within a cube of side-length 
3rvir are included, while the temperature of star forming 
gas on the effective equation of state (cEOS) is set to a con¬ 
stant value of 3000 K. The three concentric circles denote 
{0.15, 0.5, 1.0} X Tvir. 

We have separated the eight systems into two groups 
based roughly on the morphology of their halo gas - the 
first four haloes (hO-3) are closer to equilibrium at z = 2 and 
the structure of their hot gas is roughly spherical, whereas 


the second four (h4-7) are significantly more disturbed. In 
the first case, the virial shock can be clearly seen as a sharp 
increase to higher temperature and entropy, typically at 
> 1.25rvir- However, the radius of the shock (see also lSchaH 


& Springel 2015) or its existence at all depends strongly 


on direction, which we explore further in Section These 
boundaries are also associated with a decreased inflow ve¬ 
locity - that is, radial stalling - as well as increased gas 
density. 


At z = 2 haloes at this mass scale of ~ 10^^ M© typically 
reside at the intersection of multiple large-scale filaments 
(and/or sheets) of the cosmic web. This naturally leads to 
a significant amount of filamentary inflow across the virial 
radius. We see that these inflows cover only a fraction of 
the virial sphere. They can maintain coherency to at least 
half the virial radius while maintaining a strong overden¬ 
sity with respect to the mean gas density at each radius. 
They experience significant heating, typically at ~0.5rvir, 
from < 10^'^ K to > 10® K, reaching the peak temperature 
of the mean Tgas(r) in the inner halo. Their entropy can 
remain lower than that of the hot halo gas at these radii, 
^ IqS.0-8.5 K cm^ compared to ~10^Kcm^. Although the 
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Figure 3. Mass-weighted projections of gas density, temperature, entropy, and radial velocity for the first four simulated haloes at z = 2. 
In each case, all gas cells within a cube of side-length 3rvir are included, and distributed using the standard cubic spline kernel with 
h = 2.5rceii in orthographic projection. The white circles denote {0.15,0.5,1.0} rvir- Gas density is normalized to the critical baryon 
density at z = 2. The temperature of star forming gas on the eEOS is set to a constant value of 3000 K. 


mean radial velocity inside the halo is generally near equi¬ 
librium/zero (denoted by white in the colormap), we see that 
in a given direction this is rarely the case, finding instead 
either non-negligible inflow or outflow velocities. 

Two haloes, h3 and h7, have a large number of small, 
prominent, inflowing streams in the radial range 0.5rvir ^ 
r < 1.25rvir (see density maps). They originate from the 


lower left direction in the case of h3, and from the top 
for h7. These filaments are not associated with large-scale 
structure, and have characteristics distinct from the much 
larger filaments which are associated with features at larger 
radii. In particular, although they obtain similarly high in¬ 
flow velocity, their entropy and temperature are above that 
of the IGM, and their overdensity with respect to the radial 
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Figure 4. As in the previous figure, mass-weighted projections of gas density, temperature, entropy, and radial velocity at z = 2, but 
for the other four simulated haloes. 


mean is not as significant. Some of the other six haloes ex¬ 
hibit similar features at various points in the past, but they 
have disappeared by z = 2. In general, we observe that these 
features form between one and two times the virial radius 
semi-spontaneously, as in an instability, commonly triggered 
by perturbations either from substructure debris or the in¬ 
tersection of sheet-like structure in the cosmic web. They 
are marginally evident at LIO and absent at L9. We reserve 


an in-depth analysis of the formation of these features and 
whether or not their growth corresponds to a physical insta¬ 
bility in the gas for future work. 

We visualize two additional scales, the inner halo struc¬ 
ture as well as the large-scale context, of the eight simulated 
haloes at z = 2. Figure [^includes haloes 0-3 while Figure 
includes haloes 4-7. The left two columns show the projected 
density of the gas and stellar components, the white circles 
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Figure 5. Mass-weighted projections of gas and stellar density temperature at small scales (left two columns), and projected gas density 
and dark matter density at large scales (right two columns). The first four simulated haloes are shown at z = 2. All gas cells or particles 
within a cube of side-length l.Orvir (small scales) or 5.5rvir (large scales) are included, and distributed using the standard cubic spline 
kernel with h = 2.5rceii (gas) or h = 7 * 32 ,ngb (fhe radius of the sphere containing the 32 nearest neighbours of the same particle type, for 
stars and dm) in orthographic projection. The white circles denote {0.05,0.15} rvir (left two columns) and {l,2}rvir (right two columns). 
Densities are normalized to the critical (baryon) density at z = 2. 


indicating 0.05 rvir and 0.15 rvir - The right two columns show 
the projected density of the gas and dark matter compo¬ 
nents, the white circles here indicating l.Orvir and 2.0rvir- 
In the centre of the halo, r/rvir ^ 0.3, the gas com¬ 
ponent forms a complex morphology with multiple orbit¬ 


ing substructures, commonly associated with tails and gas 
streams from ram-pressure stripping, as well as tidally- 
induced bridges and spiral patterns. Rapidly varying veloc¬ 
ity fields modulate inflowing as well as outflowing material. 
At r/rvir ^ 0.05 a large gas disc is a common feature - 
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Figure 6. As in the previous figure, projections of gas and stellar density on small scales (left two columns), and gas density and dark 
matter density on large scales (right two columns) at z = 2, but for the other four simulated haloes. 


the final stage in the transformation from large-scale struc¬ 
ture and the associated tidal torques to the formation of a 
galaxy. Here the coherent angular momentum of the halo 
as a whole must transition to the angular momentum of the 
forming disc through the complicated dynamics of this inner 
halo ‘messy’ region ( [Danovich et al.|2014 ). In the small-scale 
views of h2 we see tidal features in the collisionless stellar 
distribution (extending to the right from the centre), with¬ 


out any corresponding structure in the gas. In h3, we see a 
massive tidal tail in gas extending out to ~ 0.5rvir (towards 
upper right of panel), with no corresponding structure in the 
stars. A recent post-pericenter passage of an ongoing major 
merger in h4 has induced tidal features in the large com¬ 
panion (upper right of figure) which is still on an outgoing 
orbit, and generated a long tidal tail, of which the inner half 
has reversed direction and is falling back on to the central 
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galaxy. A recent merger in h5 has left a double nucleus in 
the stars, separating by roughly one kpc. In h7 we see a 
strong tidal warp in the gas and stars of the merging com¬ 
panion (to the upper left of the central), with several stellar 
shells from earlier minor mergers superimposed. In general, 
we can see that many of the narrow inflowing streams in the 
inner halo are tidal in nature, which highlights the impor¬ 
tance of characterizing the origin of gas accretion features - 
this component of the accretion rate would have been char¬ 
acterized as a ‘stripped’ component in the methodology of 


Nelson et al. (2015). 


In the inner regions of the halos which we have identi¬ 
fied as in equilibrium, the perturbations and gas structures 


generated by the dynamics of satellite galaxies (see Zavala 


et al. 2012) are superimposed on a gas background which 


is aligned at Mpc scales with large-scale gas features. For 
example, the broad inflow in hO (towards the top of the 
image panel) with a width of ~ 50 kpc at 0.5rvir is clearly 
associated with a gas overdensity extending out to at least 
3rvir- Likewise for the broad inflow in h2 (also, top of the 
image panel) which has coalesced at the intersection of two 
such cosmic web filaments. These inflows are not smooth 
features of constant density gas. They have a ‘wavy’ or rip¬ 
pled appearance, arising from gas density contrasts as high 
as a factor of ~ 10^. By visual inspection it is clear that, at 
least at sufficient resolution, such inflowing filaments are in 
fact associations of multiple smaller structures with coherent 
kinematics and origin. 

The second set of haloes (h4-7), which are in general 
more morphologically disturbed, do not show this clear as¬ 
sociation between small-scale and large-scale gas features. 
Their time evolution indicates that this is often a conse¬ 
quence of a major halo-halo merger where the two gas reser¬ 
voirs collide and disrupt any radially coherent gas kinemat¬ 
ics as the hot halo re-equilibrates. However, recent major 
mergers are not always seen, and the nearby environment 
and topology of large-scale structure also appear to play a 
significant role in the sphericity of halo gas, independent of 
assembly history. For example, we consider the number of 
major mergers with ratio 77 ^ 1/3 experienced by each halo 
while the stellar mass is greater than M© (roughly 

2 < z < 4), as well as the maximum merger ratio of any 
merger, under the same constraint. In both cases, the mass 
ratio is determined at the time the less massive progeni¬ 


tor reaches its maximum stellar mass (see Rodriguez-Gomez 
et al.|[20T5 ). We find that hO, h2, and h3 have had no ma¬ 


jor mergers exceeding this mass ratio, and indeed none with 
?7 ^ 0.2. The outlier is hi, which did have one major merger 
(?7 = 0.8). In the disturbed group, h4 has had four major 
mergers (?7max = 0.6), h6 has had one (77 = 0.65), and h7 has 
had three (77max = 0.95). The outlier is h5, which has had 
no major mergers satisfying this criterion, the largest mass 
ratio being 77 = 0.2. 


4 THE PHYSICAL STATE OF HALO GAS 

In this Section, we quantify the gas behaviours noted above, 
first by giving the spherically-averaged profiles of gas den¬ 
sity, cell size, temperature, entropy, radial velocity, and spe¬ 
cific angular momentum, and then by examining the pair¬ 
wise correlation matrix of these same gas quantities. 


4.1 Characteristic Halo Properties 

In what follows we will normalize the distance of gas cells 
from the halo centre by the virial radius, as calculated for 
each halo at each resolution. For simplicity, we will in general 
not normalize the thermodynamic properties of gas since our 
halo mass range is so narrow. Instead, to provide a sense of 
reference, we calculate representative values for gas temper¬ 
ature, entropy, density, velocity, and angular momentum for 
a halo of Mhaio — lO^^MQat z — 2. Numerical values are 
given here. First, we take the virial temperature as 


rp _ in6-3 

2kB “ 


K 


( 2 ) 


following Barkana k, Loeb| ( |2001 ), where /x ^ 0.6 for a fully 
ionized, primordial gas. We take entropy as S' = P /where 
pressure is P = (7 — 1)Tp. A virial entropy is then defined 


= ^ = (7 - - 10® ° Kcm" 


( 3 ) 


where we have taken pvir = 200(Q6/Qrn)pcrit(^), the baryon 
fraction multiplied by two hundred times the critical density 
of the universe at that redshift. Using the chosen overdensity 
criterion, the virial radius is then 


_ f GAfhalo 

“ I 100P(z)2 


1/3 


~ 100 kpc 


in physical units, while the virial velocity is 


GAfhalo 


^ 205km/s. 


( 4 ) 


( 5 ) 


Finally, for the specific angular momentum of the halo we 
take 


= \/2A^ 


^ 10^ ° kpckm/s 


(6) 


using a spin parameter of A = 0.035 ( Barnes &: Efstathiouj 
1987 Bullock et al.|200T ). These characteristic values for the 
temperature, density, entropy, radial velocity, and specific 
angular momentum of the halo gas are typically reached at 
some reasonable fraction of the virial radius, e.g. ~(0.1 — 
0.5) X Tvir- We use them below to provide a useful reference 
with which to identify halo gas in different physical states. 


4.2 Radial Gas Profiles 

In Figure we show the median radial prohles of six gas 
quantities: temperature, density, radial velocity, cell size, en¬ 
tropy, and specific angular momentum. The main panel plots 
the profile for each of the eight haloes individually (different 
line colours) at the Lll resolution level. Light colours are 
used to indicate the four ‘equilibrium’ haloes, while darker 
colours indicate the four more ‘disturbed’ systems. The sub¬ 
panels show the ratio of the radial profiles comparing the 
L9/L11 and LlO/Lll runs, where the thick black lines in¬ 
dicate the average across haloes. We have excised all grav¬ 
itationally bound substructures before constructing mean 
radial profiles, which would otherwise appear as a forest of 
small, dense, cold systems in these figures. A few observa- 
tionally and theoretically motivated scalings with radius are 
included for reference (dotted and dashed lines). 

For instance, going beyond an isothermal model for gas 
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Figure 7. The median radial profiles of six quantities: gas temperature, density, radial velocity, cell size, entropy, and angular momentum. 
Each halo is shown separately (by colour). The main panels for each quantity show the highest resolution only, while the subpanels show 
the ratio of the two lower resolution runs to Lll (the mean ratios across all eight haloes are indicated with thick black lines). Various 
scalings with radius are provided for reference (dotted and dashed lines). All gas within 2rvir is included, except that substructures are 
excised. 


temperature, we can take a gas component tracing dark mat¬ 
ter distributed according to a NEW profile. This leads to 
a scaling of Tgas oc Menc{r)/r (e.g. Makino et al. 1998), 
which is still much too shallow (dotted line). Instead for 
r/r^ir ^ 0.15 we have approximately Tgas oc Of all 

eight systems, hO exhibits the sharpest temperature drop 
defining an outer boundary of the halo at r ~ 1.2rvir, while 


in most cases there is no sharp transition at or near the 
virial radius. In all cases the temperature drops precipi¬ 
tously in the inner halo, at 5%-10% of rvir- Note that we 
have assigned a gas temperature of 3000 K to gas on the 
star-forming equation of state, which would otherwise have 
an artificially high effective temperature. The Tgas(r’) pro¬ 
files are well converged with resolution. The mean across 
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all eight haloes and averaged over 0.15 < r/rvir < 1.5 is 
L9/L11 ~ 1.15 and LlO/Lll ^ 1.08, while at the largest radii 
the lower resolution runs have mean temperatures biased 
high by up to ~ 30%. Physics will dominate resolution - dif¬ 
ferent models for galactic-scale stellar winds or AGN feed¬ 
back can substantially modify the mean gas temperature 
profile inside as well as outside rvir ( Suresh et al.||2015 ). 

A gas density scaling as Ugas cx: approximates the 

mean behaviour of the haloes well, as does a model for an 
isothermal gas in hydrostatic equilibrium ( |Suto et al. 


using B = (n 1)Bp with n = 7, Bp = 1). However, the 
other panels make it clear that the gas is not in hydrostatic 
equilibrium, nor isothermal, implying this agreement may 
be coincidental. Convergence at lower resolutions is reason¬ 
able - over the same radial range we find L9/L11 ~ 0.8 and 
L10/L11~0.9, indicating that the mean gas densities are 
lower in the lower resolution runs. This is potentially a con¬ 
sequence of better resolving satellite galaxies and their in¬ 
teractions with the central host, which fills the halo vol¬ 
ume with more high density gas cells. These are no longer 
instantaneously bound to the satellite, and so are not ex¬ 
cluded.!^ The total halo gas mass actually decreases some¬ 
what at higher resolution (~ 10% or ~5% lower for Lll as 
compared to L9 or LIO, respectively), which would other¬ 
wise modify gas densities in the opposite direction. Finally, 
the variable assembly histories lead to the large halo-to-halo 
scatter, while the sensitivity of merger states to temporal 
offsets driven by short dynamical time-scales leads to large 
scatter between resolution levels. As with temperature, how¬ 
ever, the caveat is that radial density profiles will depend in 


detail on the implemented feedback models (Hummels et al. 


2013). 


We see a similar scatter in the Hubble-corrected gas ra¬ 
dial velocities. The haloes have different velocity structures, 
which commonly feature a slowly increasing inflow veloc¬ 
ity down to ~0.75rvir, at which point there is a noticeable 
bump towards an equilibrium value of Vrad — Okm/s. Until 
reaching the central galaxy, the radial velocity profiles are 
then roughly flat. In the following discussion of Figure 
however, we show that the mean (or median) velocity pro¬ 
file is an exceedingly poor representation of the dynamics 
of halo gas. For example, in the mean the quasi-static halo 
rarely has t'rad > 0 for 0.15 < r/rvir <1.0, while in actual¬ 
ity this value is driven down by the superposition of rapidly 
inflowing and gently outflowing components. 

The gas spatial resolution, for which we again use the 
sphere-equivalent radii rceii as a proxy, becomes better with 
decreasing radius as expected. Given our constant cell mass 
refinement criterion, it scales roughly as the cube-root of the 
density scaling. That is, rceii oc = (mceii/pceii)^^^ oc 

Pceii^ ■ sizes are a factor of four 

and two larger, respectively, scaling with the mean inter-cell 
spacing as rceii oc TboxiV“n^tof 

We overplot the gas entropy prohles with the near- 
linear scaling Sgas oc r^'^ often seen in x-ray observations 
of local clusters (e.g. George et al. 2009). Despite look¬ 
ing at the significantly more massive cluster PKS 0745-191, 


^ The filamentary feature in the upper right panel of Figure 
extending from ~0.15rvir to ~0.5rvir is an example of such a 
tidal tail. 


~ 10^^ M© in the nearby universe (z ~ 0.1), the inferred 
radial characteristics from George et al. ( |2009| are in rea¬ 
sonable agreement with the haloes simulated in this work. 
A similar level of qualitative agreement is evident at z = 2 
for ~ 10^^ Mq haloes in gas temperature, entropy, and den¬ 
sity with previous cosmological simulations at L8-equivalent 
resolution and incorporating more realistic feedback physics 
( van de Voort fc Schaye|20T2 ). In both cases, entropy rises to 
near the virial radius at which point it flattens and becomes 
roughly constant. While some haloes exhibit this gradual 
plateau (h4), the radial entropy profile of others drops as 
r~^ or faster (h2) at some radius ^ rvir- There is a cor¬ 
respondence between the equilibrium state of the halo and 
its properties beyond the virial radius - the more disturbed 
systems (h4-7, darker line colours) generally have higher en¬ 
tropy, without any strong transition. The entropy profiles, 
averaged over all eight haloes, are well converged with re¬ 
spect to L9 and LIO. 

The specific angular momentum scales roughly as 
or even more accurately as jgas oc rvc{r) given the circu¬ 
lar velocity Vc{r) oc (Menc(^)/'r)^^^ for an NFW profile of 
this halo mass. There is broad uniformity among the eight 
haloes and good convergence in the lower resolution runs, 
with a small systematic bias towards lower angular momen¬ 
tum content in the halo, L9/L11 ~ 0.88 and LlO/Lll ~ 0.96, 
again driven by high velocity tidal debris becoming better 
resolved at the highest resolution. 


4.3 Distributions Beyond Radial Dependence 

To avoid missing important features by averaging over multi¬ 
ple gas populations with distinct properties residing at sim¬ 
ilar radii, we examine the two dimensional distributions of 
these same gas quantities. Instead of only considering their 
radial dependence, in Figure we show a full pairwise cor¬ 
relation matrix, every quantity plotted with respect to ev¬ 
ery other. Here we have stacked together all eight haloes at 
Lll, normalizing only the radius of each gas cell with re¬ 
spect to Tvir of its parent halo. Despite the different thermal 
and dynamical structures of the haloes we verify that the 
key features seen individually in each are preserved in this 
collective view. To explore the various relations, we have 
separately normalized each panel ‘column by columnk That 
is, for each value along the x-axis, the colour mapping of the 
corresponding one dimensional, vertical slice is independent, 
such that each slice extends from low values (white/blue) 
through intermediate values (green/yellow) to high values 
(orange/red). Therefore, there is no global indication of the 
distribution of mass within each panel (although very noisy 
regions indicate a poor sampling and therefore comparably 
little gas mass). In addition, the corresponding panels across 
the diagonal are not just transpositions of one another, but 
answer different questions. For example, the Tgas{r) panel 
in the first column indicates, for each radius, the dominant 
temperature(s) of gas at that radius, whereas the r(Tgas) 


In the inner halo, hi and h2 have sharp drops in both Tgas and 
Sgas at a larger radius than typical, ~ 0.25rvir- In these two cases, 
while the profiles are properly centred on the galaxy, the galaxy 
is not entirely centred within the halo. This leads to more cold, 
low entropy gas beyond the typical disc size. 
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Figure 8. The mass-weighted, two dimensional correlation matrix between all pairwise combinations of six quantities: gas radius, 
temperature, density, radial velocity, entropy, and angular momentum. All eight halos are included and stacked. Each radial bin is 
normalized independently, such that the colour mapping reaches its maximum intensity at each x-value, independent of the radial 
distribution of gas mass. Therefore, panels symmetrically across the diagonal from one another are not simply transpositions, but include 
additional information. Substructures are excised. 


panel in the second column indicates, for each temperature, 
the dominant radius (or radii) of gas with that temperature. 
We describe some of the more interesting features seen in 
this matrix, which have been labelled with letters (A)-(M). 

• (*) Denotes features which arise in the intergalactic 
medium and disappear if gas at r > rvir is excluded. 
In particular, all the strong constant entropy features at 
Sgas — 10^'^ Kcm^ in the second row from the top, and the 
near-discontinuities at this same value in the second column 
from the right. 

• 1^ Regions demarcated with dotted rectangles indicate 
features arising in the galactic disc or its vicinity, and dis¬ 
appear if gas at r < 0.2rvir is excluded. This includes all low 
j gas at high density, which is already above the threshold 
for star formation in the galaxy. It also includes low temper¬ 
ature gas which has a radial velocity near zero, as well as 
the thin horizontal feature of Vr^d spanning ~0di lOOkm/s 
near r = 0 corresponding to a rotating disc. Interestingly, 
the prominent temperature bridge between the hot halo and 
galaxy at intermediate densities also largely disappears, in¬ 


dicating that a majority of this cooling occurs at r < 0.2rvir 
(see also G). 

• (A) The virial shock seen as a sharp temperature jump 
from Tigm ~ 10^'^ K to ~ 10® K. The width of this transi¬ 
tion varies greatly between the eight haloes, hO being the 
narrowest at Ar/rvir — 0.1, and as seen here in the stack 
broadened to Ar/rvir — 0.5 or more. The radius where this 
transition occurs is typically between 1.25rvir and 1.5rvir- 
After shocking, the temperature slowly increases with ra¬ 
dius, following the mean radial profile of the halo, until it 
reaches ~0.25rvir where higher densities lead to accelerated 
radiative cooling, allowing gas to join the cold ISM phase. 

• (B) The radial velocity profile of inflowing gas. This 

component speeds up as it flows in from large distances 
down to ~0.5rvir at which point it largely disappears. The 
mean velocity is approximately -75km/s at 2 rvir, -150km/s 
at rvir, and -225km/s at rvir/2. For free-fall from rest at 
infinity to a point mass of 10^^ M© we would expect a speed 
of -200km/s at twice the virial radius. The actual value is 
less due the combination of the Hubble expansion and gas 
dynamics. Given this offset at the halo outskirts, the sub¬ 
sequent scaling is roughly consistent with ^rad oc as 
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expected from free-fall, at least down to ~0.5rvir, below 
which inflow no longer dominates by mass. 

• (C) The radial velocity distribution of the hot halo gas. 
These two components overlap between 0.5 < r/rvir < 1.0 as 
gas transitions from rapidly inflowing to quasi-static. It re¬ 
ally is ‘quasi’, however, since a large mass of gas has positive 
radial velocity. Gas with r'rad > 0 arises from the dynamical 
formation of the halo atmosphere and associated splashback 
motion in the baryons (e.g. Wetzel Nagai||2014 ). In this 
stacked view the radius of this transition is notably interior 
to the virial shock, ~0.75rvir as compared to ~ 1.25rvir. In¬ 
vestigating the haloes individually we conclude that this is 
radial offset is largely a misleading feature arising from the 
stacking of radially averaged velocity profiles. Instead, the 
radius of a strong jump in temperature and entropy also 
closely corresponds to a sudden decrease of inward velocity, 
as a result of the transfer of kinetic to thermal energy. 

• (D) The virial shock seen in a sharp entropy jump 
from the nearly constant value in the IGM of ~ 10^'^ Kcm^ 
to ~ 10^ K cm^ characteristic of the hot gas for haloes of 
10^^ Mq . The entropy increase occurs over a narrower radial 
range than the corresponding temperature increase, and at a 
slightly smaller radius with respect rvir- This may be indica¬ 
tive of pre-shock compressive heating, although we caution 
that this radial difference varies from halo to halo and as 
shown in the stack. For example, in our subsequent explo¬ 
ration of hO in Figure we see that there is a close cor¬ 
respondence between the radii of temperature and entropy 
jumps. After shocking, the entropy slowly declines with ra¬ 
dius, following the mean radial profile of the halo. 

• (E) The angular momentum distribution is uni-modal 
and scatters about its mean profile, with no evidence for 
multiple populations of gas having distinct jgas at any ra¬ 
dius. 

• (F) The radial velocity depends somewhat on the tem¬ 
perature of the gas, with the hottest gas populating the 
high positive velocity tail, while cold gas is inflowing. At 
all temperatures there is a continuous distribution of radial 
velocities, with an upturn at ~ 10® K, above which gas has a 
mean radial velocity consistent with zero, and below which 
the mean I’rad is always negative (see also Joung et al.|20I2 ). 
The rapidly inflowing hot gas, discussed in (G), is the high 
velocity tail of the gas distribution at these temperatures 
and so sub-dominant by mass. 

• (G) The usual ‘phase-diagram’ plotted for cosmologi¬ 
cal simulations. Gold, low density gas (ugas ^ I0“^cm“^) 
in the IGM (denoted by a star) occupies the lower left cor¬ 
ner, the tight relation with T oc indicative of adiabatic 
compression. Shock-heated gas at intermediate densities is 
the only source for T > 10®'® K gas. At higher densities 
(^gas ^ I0“® cm“®) strong cooling flows develop at small 
radii, after which gas approaches the effective temperature 
floor of ~ 10^ K until it reaches the star formation threshold 
at riH = 0.13 cm“^. 

• (H) The heating of inflow, seen here as the dominant 
temperature for gas at each radial velocity. The handful of 
strong vertical features are largely an artefact of the stack¬ 
ing. For a given halo, when we restrict this panel to include 
only gas with 0.2 < r/rvir < I-O we find a clear correla¬ 
tion between deceleration and heating, where gas transitions 
from cold at -300 km/s to hot by -150 km/s. This is decidedly 
interior to the virial shock, and outside the galaxy. 


• (J) The average radial velocity as a function of entropy, 
restricted to gas in 0.2 < r/rvir <1.0 (not shown) resolves 
nicely into two distinct populations, where material with 
Sgas ^ lO^ ^Kcm^ has a mean rrad = 0 which is roughly 
constant. On the other hand, gas with lower entropy than 
this threshold has a mean rrad — —150 km/s, with a gradual 
trend towards faster inflow with lower entropy. 

• (K) Similarly, the average radial velocity as a func¬ 
tion of angular momentum, restricted to gas in 0.2 < 
r/r-vir < I.O resolves nicely into two components. Gas with 
j < 10^'^® kpckm/s has zero mean radial velocity, whereas 
higher angular momentum gas has t'rad — — '?^vir (this is not 
shown explicitly in the figure, which includes all gas). The 
two populations overlap for I0^‘ < jg as < 10^ ®. Together 
with (J) this panel shows the distinct physical properties of 
the quasi-static and inflowing halo gas. 

• (L) Star forming gas has a temperature on the effective 
equation of state, which is set to a low, constant value to 
differentiate it from other gas. This results in the constant 
temperature band at ~I0^® K across this entire row. This 
gas is restricted to small radii, high densities, low radial 
velocities, and low angular momenta. 

• (M) Gas with high entropy, S > lO^Kcm^ has a well 
defined relation to radius, following the mean halo profile. 
Below this threshold, if we exclude material in the IGM, gas 
at any given entropy is essentially distributed throughout 
the entire halo. 

We have seen in several cases how the mean or median 
radial profiles fail to fully capture the full state and structure 
of halo gas. Glear examples are the distributions of temper¬ 
ature, entropy, and radial velocity as a function of radius. 
On the other hand, the distributions of density and angular 
momentum as a function of radius are comparatively well 
described by a median and scatter. By measuring average 
properties within radial bins we implicitly assume that the 
structure of halo gas is spherically symmetric. Even Eigure 

fails in this regard. Eor instance, the finite radial thickness 
of the temperature jump associated with the virial shock 
could either arise from (i) a spherically symmetric feature 
of that same thickness, or (ii) the superposition of many 
thin shocks spread throughout the same radial range, dif¬ 
ferent shock fronts existing at different radii depending on 
direction. Motivated by the latter possibility, we proceed to 
investigate and quantify any variation of the structure of 
halo gas as a function of angle on the sphere. 


5 ANGULAR VARIABILITY 

As a first step away from radial averages, we consider the 
properties of gas in a broad radial shell encompassing the 
regime of interaction between quasi-static material and fil¬ 
amentary inflow. In Eigure we plot the distributions of 
temperature, density, radial velocity, and entropy, for all 
gas cells within 0.5rvir < 'Tgas < I.Orvir. In particular, we 
examine how well the lower resolution runs reproduce the 
distributions of these quantities found in the highest resolu¬ 
tion (LII) runs. In this radial range, we find that the gas of 
all simulated haloes has broadly similar properties, differing 
only in the details. 

The temperature distribution (upper left) shows a 
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Figure 9. The mass-weighted PDF of instantaneous gas temperature, density, radial velocity, and entropy, restricted to the radial 
range O.Srvir < ^gas < l-Orvir- The upper panel shows all eight haloes at the highest resolution level (Lll, different colours), while the 
subpanels below show the ratio of the two lowest resolution levels to the highest, L9/L11 (thin lines) and LlO/Lll (medium lines). The 
vertical axes of the ratio subpanels are logarithmic from 0.1 to 10. For the most part, the eight haloes show similar structure and their 
lower resolution counterparts scatter about the Lll distributions (see text for details). 


broad peak centred roughly at Tvir, slowly falling off to¬ 
wards lower temperatures and with a distinct low temper¬ 
ature peak at ^ 10^'^^ K. While the fractional amount of 
gas at each temperature in the L9 and LIO runs can differ 
throughout this regime by up to a factor of two, depend¬ 
ing on halo, the mean ratio of both L9/L11 and LlO/Lll is 
consistent with unity from 10^'^ K to 10®'^ K. At the highest 
temperatures > 10®'^ K the lower resolution runs have much 
larger deviations with respect to Lll, but we attribute this 
primarily to poor sampling of the extreme tail of the distri¬ 
bution due only to a low number of available gas cells. In 
contrast, the PDF at the low temperature peak is systemat¬ 
ically lower in the lower resolution runs. At ~ 10^'^^ K where 


there are still a large number of gas cells, the mean across 
all eight haloes is L9/L11~0.5 and LlO/Lll ~ 0.8, indicat¬ 
ing that there is a smaller fraction of the total gas mass at 
these temperatures at lower numerical resolution. Note that 
10^ K is effectively the temperature floor due to cooling in 
these simulations, and the mean Tigm at redshift two is just 
above this value. 

Gas density shows an even larger variation on a halo to 
halo basis - some have a strong second peak at higher den¬ 
sities, while this feature is largely absent for the less relaxed 
haloes. The ratio with respect to the lower resolution runs 
is consistent with unity for log Ugas < —2.5cm“^. Between 
this value and the star formation threshold of ~ 10“^ ° cm“^, 
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there is a decrease, at lower resolutions, in the fractional 
amount of gas at these densities. In this range, the mean ra¬ 
tio across all eight haloes is L9/L11 ~ 0.4 and LlO/Lll ~ 0.6. 
As with temperature, we see that the density PDFs can dis¬ 
agree between resolutions by as much as a factor of ten, but 
only in the tails when the magnitude drops to low values 
of ^10 The disagreement is similar with radial velocity, 
where the mean Ln/Lll ratios are consistent with unity for 
all values away from the extremes. 

The distribution of Urad itself is comprised of two broad 
components, roughly centred about zero, the positive and 
negative peaks indicative of outflow and inflow, respectively. 
The primary driver of inter-halo variation - for example, 
that hO has smaller minimum and maximum velocities - 
appears to be differences in assembly history, particularly a 
recent merger with another massive halo, or lack thereof. We 
explore this further in the following section. Finally, gas en¬ 
tropy behaves similarly to temperature, where we find that 
low entropy gas is strongly sub-dominant by mass with re¬ 
spect to the high entropy hot halo material. The largest 
variation between haloes occurs for log ^Sgas < lO^Kcm^, 
down to the floor at ~ 10^ Kcm^. The halo-average PDF at 
L9 and LIO reveals a smaller fraction of gas at these lowest 
entropies when compared to the Lll run. 

In general, we conclude that for any given halo, large 
variations between the three resolution levels can be seen, 
with deviations up to a factor of ten in the fractional amount 
of gas at a particular temperature, density, radial velocity, 
or entropy. These deviations occur mostly in the tails of the 
distributions which are poorly sampled at the lower reso¬ 
lution levels. Furthermore, the timing differences present in 
any single comparison likely influence some of the largest 
outliers. With respect to the average behaviour across all 
eight haloes, we find that the lower resolution runs have 
less low temperature, high density, and low entropy gas. By 
mass fraction with respect to the total gas mass in this ra¬ 
dial range, the magnitude of the effect is ~ 2 1.5) for L9 

(LIO), although it is unclear if this gas population is neces¬ 
sarily converged for these haloes at Lll. 


5.1 Structure Along Radial Sightlines 


To proceed, we consider sightlines originating outwards from 
the halo centre in different directions, which will have dif¬ 
ferent radial structure. If the hot halo gas was triaxial, for 
example, we could anticipate that the radius of a strong 
virial shock would differ along its major and minor axes. 
We therefore measure the angular variability of the thermal 
and dynamical structure of the halo by casting many such 
sight lines from the centre of each halo. Along each of these 
‘radial rays’ we sample the continuous fields of gas temper¬ 
ature, radial velocity, entropy, density, and angular momen¬ 
tum in fixed steps of Ar. Using this ensemble of rays we can 
then quantify the structure of halo gas without binning in 
spherically symmetric radial shells. 

Ray directions are set using the HEALPIX discretisa¬ 
tion of the sphere ( Gorski et al.|2Q05 ) into equal area pixels, 
which implies equal angular spacing of rays at all refinement 
levels. The total number of radial rays is Way = 
corresponding to an area subtended by each ray equal to 

firay = 4-7r/12iV^fide SI, Or 6 »ray = (180^/37riV^fide)^/^ deg. 

The sampling in the radial direction is linear and controlled 


by the parameter Wad such that Ar/r^iv = 2.0/ Wad- Our 
fiducial parameters of Wad = 400, Wide = 64 result in a 
sampling of Ar ^ 0.5 kpc and AO = 1 deg. At each point, 
mean gas properties are estimated with a tophat kernel with 
adaptive size equal to the radius of the sphere enclosing the 
Wgb = 20 nearest gas neighbours 

In Figure we show a snapshot of the combined ra¬ 
dial and angular structure of a single halo (hOLll) at z = 2. 
The twelve demarcated intervals along the y-axis denote the 
twelve base pixels of the HEALPIX discretisation. Within 
each base pixel the nested ordering scheme uses a hierar¬ 
chical quad-tree to preserve adjacency, and the four sub¬ 
intervals delineate the top nodes of each such quad-tree, 
implying that structure seen in the vertical direction is spa¬ 
tially coherent within these major and minor intervals. We 
note that there is no mass weighting or, for that matter, any 
indication of the mass distribution within each of the five 
panels, since the sampling points are smoothly distributed 
throughout the volume regardless of the underlying gas cell 
distribution. 

Focusing first on the temperature structure we clearly 
see a boundary separating the cold, intergalactic medium 
from hot, virialised gas. At this halo mass scale, the tem¬ 
perature increases by roughly two orders of magnitude, from 
K to K (dark blue to yellow/orange). In some di¬ 
rections this heating is coherent and occurs at essentially 
uniform radius (for example, 9-10, marked ‘A’), demarcat¬ 
ing a clear ‘virialization boundary’. However, across all sight¬ 
lines we also observe a temperature jump of similar magni¬ 
tude anywhere from 1.0 < r/rvir < 1-5, and smaller jumps 
can occur out to twice the virial radius. Even so, from Figure 
l^we know that this halo is one of the most spherically sym¬ 
metric, with a noticeable transition in the state of gas just 
outside the virial radius. In directions where the gas temper¬ 
ature is warm and exceeds ~10^ K out to twice rvir, we find 
a correspondence to a baryonic overdensity, with respect to 
mean at that distance. This is indicative of large-scale gas 
filaments and the heating associated with their earlier col¬ 
lapse. 

Interestingly, in directions where gas coherently pene¬ 
trates to radii smaller than the virial radius and remains 
cold (e.g. B), this same correspondence to overdensities at 
larger distance remains. Qualitatively, the existence of an in¬ 
flowing gas filament arising from the cosmic web suppresses 
a strong shock at the virialization boundary. This can arise 
from previous heating from filament formation outside of the 
halo, which increases the gas temperature to an intermedi¬ 
ate state between Tigm and Tvir- There can nonetheless be 
a (smaller) temperature jump around the mean virialization 
boundary (C), although this is not always true and gas can 
gradually heat seemingly all the way to the peak of the halo 
temperature profile just exterior to the disc region (e.g. D). 
Alternatively, this shock suppression can also arise from a 
delay of strong heating to deeper within the halo, in which 


^ This is close to the mean number of natural neighbours (or 
cell faces) of the evolved Voronoi mesh in a cosmological simu¬ 
lation ( [Vogelsberger et al.||2Q12l . Therefore the spatial scale of 
this effective smoothing is approximately matched to the scale of 
the stencil used in gradient estimation for linear reconstruction 
of fluid quantities in the MUSCL-Hancock scheme of AREPO 
(Sprmgel|20T0l. 
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Figure 10. The angular variability of gas temperature, radial velocity, entropy, density, and angular momentum at a given radius. Each 
pixel along the vertical direction represents a single radial ray, which are equally spaced in angular separation. The radial variability of 
the virial shock is evident in the shifting boundary between dark blue and yellow/orange in the left panel. Penetration of low temperature 
gas to radii smaller than rvir correlates with higher inflow velocities, lower entropy, and higher density, often extending out into the IGM 
at r ^ 2rvir- In this flgure we include only a single halo, hOLll, with A^side = 16, A^rad = 100- All substructures other than the primary 
have been excised, thereby excluding satellites. 


case there is less preheating evident at large distances. The 
radius of the strongest temperature jump then decreases, 
typically to 0.25rvir < r < 0.75rvir as in much of (E). There 
are essentially no sightlines along which the gas does not 
experience heating above Tigm at some radius. Since the ve¬ 
locity field, and so the streamlines of accreting gas, are not 
purely radial, this does not preclude the possibility that gas 
could avoid heating along its actual dynamical path. Such 
a trajectory would be curved in these panels. In some cases 
(F) we can largely rule out this behaviour, since over a large 
coherent solid angle no cold gas at r < 0.5rvir is directly 
connected to cold gas at larger radii. In general, however, 
the analysis based on radial rays is limited in this respect. 

As with density, we also see a strong correlation be¬ 
tween temperature and both entropy and radial velocity. 
Strong temperature jumps just outside the virial radius are 


also evident as sharp transitions from the gradually increas¬ 
ing (negative) IGM inflow velocity, through equilibrium, to 
a generally (positive) outflow speed indicative of the viri- 
alised halo gas. In some directions (G) multiple components 
are clearly visible, where an inner halo is comprised of hot¬ 
ter, denser, and more rapidly expanding gas which has more 
recently bounced back under its own thermal pressure. Sim¬ 
ilarly, the entropy of the low-density IGM undergoes a sharp 
increase before declining towards the halo centre. On a sight¬ 
line by sight line basis the radius of the entropy jump is 
nearly equal to that of the temperature and radial velocity 
jumps. Finally, the angular momentum of the gas is clearly 
the least affected by the virialization boundary. As the ve¬ 
locity held becomes increasingly complex towards smaller 
radii, the non-radial component contributing to the angular 
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momentum does not have a clear correlation with any of the 
other local gas properties. 

5.2 Different Gas Heating Regimes 

To understand the properties and the importance of sight¬ 
lines exhibiting different radial behaviours, we would like to 
identify a few characteristic types. We separate rays into a 
small number of disjoint sets in Figure ED each group ide¬ 
ally having a distinct radial behaviour. By visual classifica¬ 
tion we identify four such groups, which together encompass 
the majority of sight lines. In particular, we split rays into 
‘distant shock’, ‘closer shock’, ‘gradual heating’, and ‘low 
entropy’ types, each of which is based on a quantitative se¬ 
lection applied across some radial range. Specifically, for the 
‘distant shock’ type we require all of 

• max(dT/dr|o. 2 <r< 2 .o) > 0.25 

• 1.5 ^max(dT/dr) ^ 2.0 

• min(Tgas|o.2<r<i.4) > 2 X lO^K 

where dT/dr denotes the derivative of gas temperature with 
respect to radius, in units of (log K)/(0.1 rvir), 'rmax(Q) indi¬ 
cates the radius where the quantity Q reaches its maximum 
value, and Q|ri<r<r 2 indicates that a quantity Q is con¬ 
strained only within the radial range between ri and r 2 . For 
the ‘closer shock’ type we require 

• max(dT/dr|o. 2 <r< 2 .o) > 0.25 

• 0.8 ^max(dT/dr) ^ 1*5 

• min(Tgas|o.2<r<o.7) > 3 X 10^ K. 

That is, both must exhibit a strong temperature jump at ei¬ 
ther large or intermediate radii, while excluding rays which 
subsequently drop to low temperature at smaller radii - an 
artefact of intersecting non-radial cold debris, as we subse¬ 
quently discuss. Of the remaining rays not meeting the prior 
two conditions, we further require rays of a ‘gradual heating’ 
type to satisfy 

• max(dT/dr|o. 2 <r< 2 .o) < 0.8 

• niin(Tgas 1 0.2<r <0.75 ) ^ 3 X 10 K. 

Finally, the ‘low entropy’ type requires of any remaining rays 

• max(5'gas|o.2<r<2.o) < 8 X lO^Kcm^. 

In Figure pr] we show all four groupings as separate panels 
(for halo hOL9). The temperature is plotted as a function of 
radius, while colour indicates gas entropy. A single prototyp¬ 
ical ray is included as a thick line, while all rays belonging 
to that type are shown underneath as thin lines. The mean 
radial velocity profile of all rays of that type, locally aver¬ 
aged in radius, is denoted by the series of numbers shown in 
each panel, in units of km/s (negative denoting inflow). The 
top inset above each panel plots the temperature deriva¬ 
tive dT/dr, positive denoting increasing temperature with 
decreasing radius, for the prototypical ray (thick) and five 
other random sightlines of that type (thin). 

We briefly describe the behaviour of each type. Sight¬ 
lines experiencing a ‘distant shock’ undergo a jump in tem¬ 
perature from ~ 10^ K to ~ 10® K over a radial range of 
~0.1rvir (although, as in the example, multiple jumps can 
exist and be spread over a larger radial range). At this same 
radius the entropy also increases by approximately two or¬ 
ders of magnitude, from ~ 10^ ® K cm^ to ~ 10^ ® K cm^. The 


Table 3. The percentage of radial rays of each of the four types: 
‘distant shock’, ‘closer shock’, ‘gradual heating’, and ‘low en¬ 
tropy’. The mean fraction across all eight haloes is calculated 
separately for each resolution level. Since rays cover equal angu¬ 
lar area, each fraction corresponds to the geometrical percentage 
of the sphere occupied by sightlines satisfying each criterion. The 
errors represent the standard deviation among the eight haloes. 


type 

L9 

LIO 

Lll 

distant shock 

23.6 ± 4.8 

17.3 ± 4.9 

12.5 ± 6.1 

closer shock 

22.5 ± 15.9 

27.2 ± 14.0 

30.4 ± 15.2 

gradual heating 

28.0 ± 15.2 

14.5 ± 7.9 

6.2 ± 3.9 

low entropy 

7.1 ± 2.7 

7.5 ± 3.8 

5.4 ± 2.9 


rapid inflow velocity decreases, reaching its maximum (pos¬ 
itive) value near the radius of maximum temperature. To¬ 
wards smaller radii, the ray temperature then follows the 
steadily increasing mean temperature prohle of the halo, 
until reaching the centre. The ‘closer shock’ sight lines have 
the same behaviour - the radial distinction between close 
and distant shocks is arbitrary. In either case, the temper¬ 
ature profile can be non-monotonic if the shock occurs at 
sufficiently large radius, such that the shock temperature is 
greater than the local mean Tgas(7’), and the gas can subse¬ 
quently cool. Heating can be associated with one or multi¬ 
ple shocks (e.g. two, in the case of the prototypical example 
shown in the upper left panel). The maximum derivatives of 
temperature and entropy are ~ 2.0 (in log) per 0.1 rvir, al¬ 
though the typical maximum along any given ray is roughly 
half as large. 

The ‘gradual heating’ sightlines generally reach the 
same maximum temperature (in the inner halo) and en¬ 
tropy (near the virialization boundary), but do so without 
any sudden jumps. Here the inflow velocity is roughly con¬ 
stant and always negative, never approaching a quasi-static 
state, increasing towards smaller radii and with a mean of 
approximately -75 km/s. The ‘low entropy’ rays are selected 
to have a maximum entropy of less than 8 x 10®Kcm^ al- 
though the exact value is arbitrary. In this case, the inflow 
velocity is always strongly negative, peaking in the inner 
halo, with a mean of approximately -135 km/s. The vast 
majority of these rays still reach high temperature, but at 
<0.25rvir and possibly not until the disk-halo interface. At 
r > 0.5rvir we see that ‘low entropy’ sightlines also have sys¬ 
tematically lower temperatures and higher inflow velocities 
than the other three types. 

In Table [3] we include the fraction of each of these 
four ray types, calculated as the mean over all eight haloes. 
Although the balance between the two shock types shifts 
with resolution, the total ray fraction experiencing a strong 
shock remains at ^ 45%, with the temperature jumps mov¬ 
ing somewhat inwards. The gradual heating fraction drops 
sharply at higher resolution levels, either because fluctua¬ 
tions become better resolved, or because the adopted cut 
on the temperature derivative is slightly too restrictive. The 
low entropy ray fraction remains fairly converged at ~ 6%. 

We have verified that these values, and all the other 
quantitative results related to the radial rays, are well con¬ 
verged with the numerical parameters Wad = 400, Wide = 
64. Given our selections, approximately ~ 15% of the sight- 
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Figure 11. The temperature profiles of individual radial rays (four main panels), with colour indicating gas entropy. Each of the four 
panels includes a disjoint subset of the entire ray set, where the selection for each was chosen by visual inspection in order to find four 
types of sightlines, each with similar radial properties, which together cover the majority of behaviours. In each panel, all rays are shown 
as thin lines, while a single prototypical example is shown as a thick coloured line. The mean radial velocity profile for each ray type 
is indicated by the series of numbers (in units of km/s, negative denoting inflow). The smaller top insets show the derivative of gas 
temperature with respect to radius, in units of (log K)/(0.1 rvir)- The percentage in each panel indicates the fraction of all rays satisfying 
the criteria, for this particular halo. Given our selections, some sightlines are excluded from all four types. The vast majority of these 
rays either fall just outside a selection, or experience a significant temperature dip within the halo, which is a signature of satellite debris 
or intersection with non-radial filamentary inflows. Two such examples are shown as thin black lines in the upper two panels. In this 
figure we include only one halo, hOL9, with N^ide = 8 cind A^rad = 100 for visual clarity. 


lines remain excluded from all four types at L9, although 
this increases to ~45% at Lll. The vast majority of these 
excluded rays either (i) fall just outside a selection, or (ii) 
experience a significant temperature dip within the halo. By 
definition they do not satisfy the ‘gradual heating’ or ‘low 
entropy’ conditions. Two such examples are shown as thin 
black lines in the upper two panels. We find that the narrow 
dips are typically a signature of intersecting satellite debris 
or satellite outskirts, which were not calculated as locally 
gravitationally bound and so were not excised. Gas cells of 
this type fill a larger fraction of the halo volume at higher 
resolution, leading to the increased percentage of rays ex¬ 
cluded from all four types. Broader dips are typically inter¬ 
sections with large filamentary inflows which are not aligned 
in the radial direction. This exposes the main caveat of the 
above analysis - namely, that gas inflow with a tangential 


velocity component need not evolve according to the temper¬ 
ature, entropy, or velocity structure of any particular radial 
ray. 


5.3 Quantifying The Asphericity of Temperature 

In order to quantify the instantaneous halo temperature 
structure, we calculate differential histograms of radii satis¬ 
fying a given criterion in Figure(top panel). First, solid 
lines show the distribution of the maximum radius, for each 
radial ray, where the temperature exceeds some fraction of 
Tvir (different colours). To avoid washing out features by 
stacking, we include only one halo (hOLll). We see that for 
the three lowest temperature cuts, (0.1 — 0.5) x Tvir, the 
^max values all show a strong peak at ~ 1.2 rvir, indicative 
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Figure 12. (Top) The histogram of radii satisfying three differ¬ 
ent types of criteria (different linestyles), computed separately for 
each radial ray. Since each subtends equal angle, this histogram is 
equal to the fractional solid angle of the whole sphere covered by 
rays satisfying each criteria, as a (differential) function of radius. 
The criteria are: (i) the maximum radius at which the temper¬ 
ature along the ray exceeds some fraction of Tvir, (ii) likewise, 
but the minimum radius, and (iii) the radius at which each ray 
reaches its maximum temperature. For the first two conditions, 
we consider six different temperature thresholds (from blue to 
red). In this figure we include only one halo at the high resolu¬ 
tion level, hOLll, with A^side = 64 and A^rad = 400. (Bottom) As 
in the upper panel, but here as a cumulative function of radius. 


of the virialization boundary. The distributions are broad, 
however. In the case of T/Tvir > 0.5 the spread of maxi¬ 
mum radii is roughly symmetric, extending from 0.5rvir out 
to l.Grvir- That is, although the majority of sight lines first 
exceed 0.5Tvir at a well-defined surface sitting at 1.2rvir, a 
non-negligible fraction first exceed this temperature already 
by l.Grvir, and an even larger fraction cross this threshold 
inside the virial radius, between 0.5 — 0.75rvir. 

The three higher temperature thresholds (green, or¬ 
ange, red) reflect the increasing mean radial temperature 
profile of the halo. The virial temperature is typically 
reached at ~ 0.5 Tvir, while gas reaches twice Tvir just prior 
to the galaxy at ^0.1 Tvir- The dashed lines show the min¬ 
imum radius where each temperature threshold is reached. 
They are uniformly peaked at the halo centre. Finally, 
the dotted line plots the distribution of the radius where 
each ray reaches its maximum temperature, which are also 
peaked in the halo centre. One of our principal goals is to 


identify gas shocking at the halo-IGM transition. We see 
that the peak of the radial distribution, for rays satisfying 
r’max(T/Tvir > {0.5, 0.2}), measures the most dominant ra¬ 
dius for this transition. With respect to halo to halo vari¬ 
ability, we find that the shapes and overall widths of each 
r’max distribution can vary substantially. For instance, hO 
has a strong peak which is both the narrowest in radius 
and the most covering in angular fraction. Yet, it also has 
a prominent secondary peak of roughly half the geometrical 
importance at larger radius, Ar ~ 0.2r/rvir further from the 
halo centre. Even more extreme, h3 has two peaks of equal 
strength spaced Ar ~ 0.5r/rvir apart. Multiple radial peaks 
are a common occurrence. 


The widths of the distributions also encode a measure 
of the spherical symmetry of any strong heating. Here it is 
useful to look at the cumulative histograms of radii satisfy¬ 
ing these same temperature criteria, as shown in the bottom 
panel of Figure The slope of this CDF indicates angular 
uniformity - presence of a Heaviside step function would in¬ 
dicate perfect spherical symmetry, while slopes approaching 
zero would indicate that the temperature transition takes 
place, depending on sightline, over widely disparate radii. 
The fact that several temperature thresholds plateau below 
unity indicate that only ~85% of rays reach T > 2.0Tvir 
(red line) and ~ 95% of rays reach T > 1.5Tvir (orange line), 
while all other temperature thresholds are exceeded along 
all sight lines. Any value for the cumulative solid angle can 
provide a good way to quantify the location of the mean viri¬ 
alization boundary. For instance, the radius by which 27r sr 
are covered by rays with maximum temperature exceeding 
half of Tvir is ~ l.lrvir- 


We conclude with a visual impression of the three- 
dimensional temperature structure in and around a single 
halo (hOLll) in Figure 13 A ray-traced volume rendering 
highlights five iso-temperature surfaces using the same blue- 
red colortable and physical bounds of 4.4 < log Tgas [K] < 
6.3 as in the temperature projections of FigureIt is sam¬ 
pled with a transfer function comprised of four narrow alpha- 
channel gaussians at log Tgas [K] = {4.5, 5.0, 5.5,6.0} with 
widths of cr = {6,6,3,3} x 10“^ in log K. The scale of the im¬ 
age is 5.5 Tvir in both height and depth, as in the zoomed out 
panels of Figure The dark blue surface straddles the mean 
Tigm at ~ 2rvir and reveals a morphology at overdensities 
just below cosmic web filaments. The light blue surface at 
10^ K clearly outlines the structure of Mpc-scale, large-scale 
filaments, while the yellow surface at 10^'^ K is restricted 
to the interior of warmer filaments. Both of these two iso¬ 
temperatures surfaces are closed. Further, they enclose the 
virialised halo itself - we can see how the filamentary shape 
of the 10^ K surface at r > 2rvir transitions smoothly into a 
quasi-spherical shape by opening up, wrapping around and 
encompassing the higher temperature shells. The mean ra¬ 
dius of the 10^'^ K surface is ^l.Orvir, while that of the 
10® K surface is ~0.5Tvir- The intersection of radial sight¬ 
lines cast outwards in different directions with this complex 
temperature structure in and around the 10^^ M© halo gives 
rise to the four distinct ray types described earlier. 
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Figure 13. A volume rendering of the three-dimensional temperature structure around a single halo (hOLll) at z = 2. We use an 
orthographic projection with depth and height equal to 5.5rvir (—620kpc), as in the zoomed out panels of Figure]^ The colortable and 
4-gaussian transfer function are shown. A front-to-back ray tracing method is used, with no scattering and no absorption. Along each 
ray we sample with a constant step-size of 0.25 kpc. The temperature at each sample point is calculated using the standard cubic-spline 
SPH kernel interpolant with adaptive smoothing length h over the N = 200 nearest neighbours (the representation of the gas in terms 
of a Voronoi tessellation is not used). The white dot marks the halo centre, and the white circles denote one and two times the virial 
radius. 


6 DISCUSSION 


6.1 Observational Points 


We have deferred a direct comparison against observations of 
the gas content of haloes to future work. Such a comparison 
deserves a careful treatment of the various steps required 
to make robust synthetic (or ‘mock’) observations from the 
simulations, for both hydrogen (e.g. [Faucher-Giguere et aT7 
20111 IFaucher-Giguere et al. 112014] |Fumagalli et al.[ 2014| 
Bird et al.|2014 ) as well as metal signatures (e.g. Shen et al. 


2013 [Hummels et al.|[2013| [Ford et aL]|2014| jSuresh et al. 


2015). As motivation and in preparation for the next work 


in this series, there are a number of insightful observations 
of the gas in and around haloes with masses of ~10^^ M© at 
zr^2 which prompt mention here. 

Observations have studied the incidence and kinematics 
of HI absorption around haloes of this mass (e.g. the Lyman- 
break galaxies of Rudie et al.||2012 2013). They found, for 
instance, a covering fraction for Lyman-limit systems within 
the virial radius of ~ 30%. The hydrogen and metal content 
around lower mass damped Lyman-alpha systems (DLAs, 
~ 10^^ Mq haloes) at this re dshift is roughly con sistent with 
their 10^^ M© counterparts (Rubin et al. 2014), where the 


covering fractions within 100 kpc of strong Si II absorption 
was found to be ~ 20%, versus ~ 60% for strong C IV. 

Around more massive ‘quasar hosts’ of ~ 10^^'^ M© the 


et al.||2013 ). On the scale of the virial radius in such haloes 
the covering fractions of cold metal ions (e.g. C II and C IV 
Prochaska et al.|2014 ) are even higher and may approach 


unity. The presence of cold gas at large radii (Turner et al. 


2014b) is a result which, in general, current hydrodynam- 


ical simulations have significant difficulty reproducing. On 
the other hand, hot gas as often probed by O VI at these 


redshifts (Turner et al. 2014a) is more readily reproduced 
given sufficiently energetic outflows ( Ford et al.|2013 Suresh 
et al.|20l5 ). In some cases there is strong evidence that the 
observed absorption results from high velocity galactic out¬ 
flows ( Crighton et al.||2015 ). Other observations imply that 
slower inflowing structures with a high degree of kinematic 
coherence are more likely ( Rubin et al.|2014| . In general, it is 
unclear if the cold gas seen in the haloes surrounding galax¬ 
ies is always the consequence of one of (i) feedback-driven 
outflows (e.g. Marinacci et al. 2014b| or (ii) cosmological 
inflow, or if this depends on the halo mass under investiga¬ 
tion. Either way, it is currently difficult, from the theoretical 
perspective, to explain how this cold gas is generated, dis¬ 
tributed, and either maintained or replenished in the halo. 


6.2 The Resolution Issue 

Observations have also placed constraints on the physical 
size of the systems which give rise to metal absorption in 


covering fractions of neutral hydrogen are larger (Prochaska gaseous haloes. Simcoe et al. (2006) estimated absorbers 
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sizes around z ^ 2.3 galaxies ranging from the sub-parsec 
to the kiloparsec. This is consistent with the transverse 
scales of C IV absorbers derived in Rauch et al.| (2001) 
from multi-sightline analysis of lensed QSOs at sub-kpc 



zc^2.b galaxy and derive a (< 100 — 500) pc size constraint, 
which they argue is a ubiquitous size-scale for all low-ion 
halo gas in the circumgalactic environments of both low 
and high redshift galaxies. Fieri et al. (2014) conclude that 
a typical low-ionization absorber associated with the CGM 
of strong Lyman-alpha forest systems at z ~ 2.5 is ~30pc 
in size. Furthermore, the gas giving rise to various absorp¬ 
tion lines may also have unexpected geometries with little 
resemblance to spherical clouds (Churchill et al. poTIl i. 

Obtaining ~ 100 pc hydrodynamic resolution at the 
virial radius by simply running higher resolution simulations 
of isolated haloes will be prohibitively expensive. Following 
the observed scaling of the gas resolution in the halo, we 
would obtain a mean rceii — 0.1 kpc physical at z = 2 with a 
LI 4-class simulation. That is, gas mass resolution approach¬ 
ing 30 Mq , requiring > 10^° gas cells in a 10^^ M© halo. This 
is not a realistic goal for the near future. Even given such 
a simulation, the minimum size of resolved clouds would be 
a few times larger than the size of individual resolution ele¬ 
ments. 

On the other hand, numerical methods may enable 
novel ways to focus resolution within the halo regime. If we 
can avoid clustering the computational effort in the densest 
regions of space - within galaxies themselves - much higher 
spatial resolutions can be achieved in lower density media. 
In AMR simulations, the refinement criterion can be chosen 
as desired, on the gradient of the neutral hydrogen fraction 
for instance (Rosdahl & Blaizot 2012), or with a progres¬ 
sively refined uniform Eulerian grid which is volume filling 
( Miniati|2014 ). In our moving mesh simulations, we can en¬ 
force a target gas cell mass criterion of any type through 
adaptive refinement and de-refinement of Voronoi cells. In 
future work we will explore this possibility as a solution to 
the resolution issue, in addition to quantifying the rates and 
modes of accretion using the tracer particle information. 


7 CONCLUSIONS 

In this first paper of the series we present a suite of high- 
resolution cosmological hydrodynamic zoom-in simulations 
targeted at understanding the properties of halo gas and 
cosmological accretion in ~ 10^^ M© haloes at z = 2. Using 
the moving-mesh code AREPO we simulate each of eight 
haloes at three levels of increasing resolution, reaching a 
mean baryonic mass of ~10,000 solar masses. We study the 
thermal and dynamical state of halo gas within 2 r^ir and in 
particular quantify density, temperature, entropy, angular 
momentum, and radial velocity in terms of their mean radial 
profiles and, more interestingly, in terms of their variation 
along radial sight lines in different directions. 

In this mass regime, haloes typically reside at the inter¬ 
section of one or more ‘cosmic web’ filaments arising from 


large-scale structure. The result is a significant amount of 
filamentary inflow across the virial radius. 

• Investigating the interaction of this inflow and the 
quasi-static hot halo atmosphere, we find that stream mor¬ 
phologies become continuously more complex with better 
numerical resolution. In general, single coherent flows tend 
to resolve into multiple, narrower streams while producing 
density and temperature structure at smaller spatial scales. 
However, we point out that even at our highest resolution - 
comparable to the best zoom-in simulations that currently 
exist for haloes of this mass - the gas-dynamics in the cir¬ 
cumgalactic regime are poorly resolved in comparison to 
within galaxies themselves. 

• In general, the mean radial profiles of gas properties are 
well captured at our lowest resolution - equivalent to that 
currently available in full cosmological volume simulations. 
However, gas which is stripped from infalling satellites, or 
which forms as the result of a strong tidal interaction, is less 
well converged, and populates the halo with additional high 
density gas in our highest resolution runs. 

• Although average radial profiles are well converged, 
their direct interpretation can be misleading. Within the 
halo itself, 0.2 < r/rvir < 1-0, there are clearly multiple 
gas components overlapping in radius and corresponding 
to quasi-static versus inflowing material. Eor example, al¬ 
though the mean radial velocity within the halo is near zero 
(equilibrium), this is not true for the majority of gas, which 
is actually either inflowing or outflowing. 

• Examining the thermal and dynamical state of gas 
within 2rvir, we clearly identify the existence of a strong 
virial shock. This ‘virialization boundary’ typically resides 
at (1.25 — 1.5) X Tvir and is evident as a sharp increase in 
gas temperature, entropy, density, and a decrease in inwards 
radial velocity. Collimated inflows which remain cold past 
the virial radius experience significant heating, typically at 
~ O.Srvir, from < 10^'® K to > 10® K. 

• Although the mean radius of a strong virial shock may 
be ^ 1.25 Tvir, our sightline analysis allows us to identify, in 
each halo, many radially-sharp virial shocks, each at some 
radius between l.Orvir and 1.5rvir, depending on angle on 
the virial sphere. 

• Investigating the process of gas virialization, we identify 
different mechanisms responsible for the heating of gas in the 
circumgalactic regime. In addition to a single strong virial 
shock at ^Tvir, we find that gas can shock at much larger 
radii, <2rvir, particularly in systems where the hot halo is 
far from equilibrium and can extend to these distances due 
to the dynamical response of a recent major merger. We also 
see that gas can heat gradually, with temperature increasing 
slowly from its characteristic IGM value until reaching the 
mean halo temperature profile. Einally, we identify the ex¬ 
istence of radial sightlines along which gas entropy remains 
always below the level characteristic of the hot halo. Even in 
this case, however, gas cannot avoid heating to ~ Tvir in the 
inner halo, ~ 0.2rvir, just prior to cooling onto the centrally 
forming galaxy. Both of these two cases - gradual heating 
and persistent low entropy - are sub-dominant, accounting 
geometrically for only ~ 15% of the total 47r, with a strong 
shock covering the other ~ 85%. 

• Einally, we conclude by assessing the ashpericity of halo 
gas and by measuring the location and width of the virial- 
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ization boundary as demarcated by temperature. We find 
that the distributions of maximum radius rmax, where the 
temperature exceeds e.g. (0.1 — 0.5) x Tvir along each ra¬ 
dial sightline, can be quite broad. In general, for the less 
disturbed systems, the majority of radial sight lines first sur¬ 
passes 0.5 Tvir at a well-defined surface sitting at or just 
beyond rvir- However, a non-negligible number already ex¬ 
ceeds this temperature at larger distance, and an even larger 
fraction crosses this threshold only inside the virial radius, 
between 0.5 — 0.75rvir. We propose that this geometrically- 
motivated analysis can be used to quantify the structure of 
gas heating as a result of the halo-IGM transition. 

In order to serve as a benchmark for the realization of 
these eight haloes, we have intentionally not included ex¬ 
isting models for energetic feedback due to star formation 
(in the form of a kinetic galactic-scale wind generation) and 
AGN (in the form of quasar and radio mode thermal en¬ 
ergy input plus local radiative effects). We have deferred 
any direct comparisons with observations of the gas content 
of haloes. Both points remain directions for future work, 
particularly in the sense of how simulations without galactic 
winds can reproduce observed hydrogen and metal signa¬ 
tures. The answer may depend sensitively on details of the 
numerics, including hydrodynamic resolution, particularly if 
the spatial scale of dense gas structures in galactic haloes is 
as small as commonly estimated. 
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